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EXECUTIVE  SUMMARY 


This  paper  serves  to  elucidate  the  mathematical  steps  needed  to  apply  near-field  imaging  over  a 
significant  (angular)  area  of  regard,  with  specific  application  to  the  receiving  sites  of  the  Air  Force 
Space  Surveillence  System  (AFSSS,  “Space  Fence”).  Cornwell  has  developed  the  mathematical 
formalism  to  merge  near-field  imaging  with  imaging  over  a  wide  field  of  view.  The  traditional 
far-held,  narrow  held-of-view  imaging  developed  in  radio  astronomy  can  be  extended  to  wide  fields 
of  view  by  taking  into  account  the  three-dimensional  nature  of  the  receiving  array  (the  so-called 
“re  term”),  leading  to  a  method  called  “re- projection.”  This  re-projection  method  has  the  effect 
of  retaining  a  Fresnel  term  in  the  imaging  process,  and  Cornwell  sketched  how  this  re-projection 
might  also  be  applicable  to  near-held  imaging.  For  the  Space  Fence  receiving  sites,  the  near  held  is 
located  an  altitude  of  approximately  1000  km.  I  illustrate  how  a  near-held,  re-projection  imaging 
could  be  implemented  for  determining  the  direction  cosines  of  radar  reflections  from  a  system  like 
the  Space  Fence. 
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ON  NEAR-FIELD  LU-PROJECTION  FOR  RADIO  INTERFEROMETRIC 

IMAGING 


1.  BACKGROUND 

In  traditional  radio  astronomical  interferometric  observations  of  celestial  sources  outside  the 
solar  system,  imaging  can  be  performed  under  the  far-held  approximation  (Thompson,  Moran,  & 
Swenson,  1986,  Chapter  3).  For  objects  within  the  solar  system,  however,  near  field  imaging  may 
be  required.  Specific  examples  from  radio  astronomy  include  Very  Long  Baseline  Observations 
of  radar  returns  from  near-Earth  objects  (Black  et  al.,  2005).  Although  various  techniques  were 
developed  to  handle  near-field  imaging  have  been  developed,  we  believe  that  Carter  (1988)  was 
the  first  to  place  these  on  a  firm  mathematical  footing.  In  particular,  Carter  (1988)  distinguished 
between  the  near-held  of  the  array  and  the  near-held  of  the  source. 

While  radio  astronomical  imaging  of  objects  in  the  solar  system  may  require  taking  near- 
held  effects  into  account,  the  density  of  objects  in  the  sky  is  sufficiently  low  that  searches  at 
radio  wavelengths  (with  the  sensitivities  provided  by  current  instrumentation)  are  not  practical. 
Consequently,  only  a  small  region  of  the  sky  surrounding  the  target  object  needs  to  be  imaged.  In 
contrast,  one  of  the  operational  modes  for  a  space  situational  awareness  array,  such  as  the  Air  Force 
Space  Surveillance  System,1  (henceforth  termed  the  “Space  Fence”)  involves  searching  for  radar 
returns  from  spacecraft.  A  large  region  of  sky  (e.g.,  >  90°)  may  need  to  be  searched  efficiently. 

More  recently,  Cornwell  (2004a)  has  extended  the  work  by  Carter  (1988),  making  use  of 
the  w- projection  methodology  (Cornwell  et  al.,  2008).  Briefly,  the  measured  visibilities  or  cross¬ 
correlations  from  a  radio  astronomical  interferometer  are  intrinsically  3-dinrensional  while  the  de¬ 
sired  sky  brightness  is  a  2-dinrensional  quantity,  i.e.,  a  function  of  angles  or  direction  cosines.  In 
standard  (far-held)  radio  astronomical  processing,  particularly  for  narrow  helds  of  view,  one  can 
approximate  the  measured  visibilities  by  a  2-dinrensional  quantity  and  then  use  a  Fourier  trans¬ 
form  relation  to  obtain  the  sky  brightness.  This  approximation  is  equivalent  to  assuming  that  the 
interferometer  is  essentially  planar.  The  ic-projection  relaxes  these  approximations  by  making  use 
of  higher  order  terms  that  are  neglected  during  the  standard  radio  astronomy  processing.  Essen¬ 
tially,  the  difference  between  standard  far-held  and  tc-projection  is  that  the  standard  method  uses 
only  the  Franhaufer  term  while  the  rc-projection  retains  a  Fresnel  term.  Using  a  Fresnel  term  for 
re-projection  also  implies  that  it  might  have  applicability  to  near-held  imaging  over  large  regions 
of  sky. 

This  document  serves  to  elucidate  various  assumptions  or  steps  that  Cornwell  (2004a)  leaves 
unstated,  with  the  particular  aim  of  applying  re-projection  to  the  Space  Fence.  Throughout  we 
use  the  same  coordinate  system  (Figure  1)  as  Cornwell  (2004a):  The  positions  of  receptors  or 
antennas  are  denoted  by  r,;  (i.e.,  numerical  subscripts),  the  positions  of  emitters  are  denoted  by  rx 
(i.e.,  alphabetical  subscripts),  and  the  direction  from  the  center  of  the  array  and  the  center  of  the 
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1This  system  was  formerly  known  as  the  Navy  Space  Surveillance  (NAVSPASUR)  system. 
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Fig.  1 — Geometry  for  u>-projection.  Receptors  are  located  at  V\  and  r2,  an  emitter  is  located  at  r0,  and  the 
distance  from  the  center  of  the  array  to  the  center  of  the  field  of  view  is  D.  All  vectors,  ri,  r2,  and  r0,  are 
three-dimensional  (i.e.,  their  2  components  are  not  assumed  to  be  0). 


(angular)  area  of  regard  is  h.  The  distance  from  the  array  of  receptors  to  the  mid-point  of  the  area 
of  regard  is  D.  We  shall  make  the  restriction  that  Zi  <  D  and  zx  <  D. 

2.  PHYSICAL  MOTIVATION  AND  RADIO  ASTRONOMY  BACKGROUND 

We  motivate  the  remainder  of  this  document  (as  well  as  the  discussions  in  Carter  1988,  Cornwell 
2004a,  and  Cornwell  2004b)  by  considering  a  series  of  increasingly  complicated  scenarios  (viz. 
Figure  1). 


1.  The  emitter  is  in  the  far- field,  and  area  of  regard  is  small.  Thus,  D  00,  small-angle  approxi¬ 
mations  can  be  used  (with  respect  to  ra),  and  the  wavefronts  from  the  emitter  can  be  assumed 
to  be  planar  across  the  entire  receiving  array.  This  case  is  the  traditional  radio  astronomical 
interferometric  imaging  case  (Thompson,  Moran,  &  Swenson,  1986).  A  fundamental  theorem 
(the  van  Cittert-Zernike  theorem)  for  this  case  is  that  there  is  a  2-dimensional  Fourier  relation 
between  the  brightness  or  intensity  distribution  on  the  sky  and  the  second  moment  of  the 
electric  field  across  the  receiving  array  (the  visibility  function) .  A  key  assumption  of  the  van 
Cittert-Zernike  theorem  is  that  the  receiving  array  can  be  assumed  to  be  planar. 
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2.  The  emitter (s)  is  in  the  far- field,  but  the  area  of  regard  is  not  small.  In  this  case,  D  — >  oo 
and  the  wavefronts  from  any  particular  point  in  the  area  of  regard  are  planar  across  the 
receiving  array.  Effectively,  the  geometry  of  the  receiving  array  as  seen  from  the  area  of 
regard  must  be  taken  into  account  or  the  non-coplanar  nature  of  the  receiving  array  can 
no  longer  be  neglected.  The  full  three-dimensional  characteristics  of  the  visibility  function 
must  considered.  The  rc-projection  methodology  was  developed  to  handle  this  case  for  radio 
astronomical  imaging  (Cornwell  et  al.,  2008);  see  also  Cornwell  &  Perley  (1992)  for  alternate 
approaches  to  this  case. 

3.  The  emitter  is  no  longer  in  the  far-held,  but  the  area  of  regard  is  small.  In  this  case,  the 
curvature  of  the  wavefronts  across  the  array  must  be  taken  into  account,  but  small- angle 
approximations  can  often  be  retained.  This  case  was  considered  by  Carter  (1988). 

4.  The  emitter(s)  is  no  longer  in  the  far  held  and  the  area  of  regard  is  not  small.  In  this  case, 
both  the  curvature  of  the  wavefronts  from  any  point  in  the  area  of  regard  and  the  non-coplanar 
nature  of  the  receiving  array  must  be  taken  into  account.  This  case  is  considered  by  Cornwell 
(2004a)  and  in  this  document. 


Finally,  we  adopt  nomenclature  and  terminology  derived  from  radio  astronomical  heritage. 
Specifically,  the  visibility  function  V  has  multiple  meanings.  From  considerations  of  the  van  Cittert- 
Zernike  theorem  and  its  extensions  to  wide-held  and  near-held  imaging,  the  visibility  function  can 
be  used  to  indicate  the  continuous  function  dehned  everywhere  in  the  3-dinrensional  volume  dehned 
by  the  receiving  array.  In  general,  the  receiving  array  consists  of  only  a  hnite  number  of  receptors, 
located  at  r*.  Thus,  the  measurements  of  the  visibility  function  are  hnite  in  number  and  occur 
only  at  set  of  locations  dehned  by  the  receptor  positions  {r*},  denoted  V({ri}).  Finally,  a  specihc 
measurement  of  the  visibility  function  is  determined  by  any  two  receptors,  denoted  V(ri,rj). 

3.  “GEOLOCATION  AS  A  SEARCH  TECHNIQUE” 


Cornwell  (2004a)  casts  geolocation  as  a  search  technique  for  hnding  the  source(s)  of  emission 
from  the  received  voltage  streams  from  an  array  of  antennas.  The  “visibility”  for  a  pair  of  antennas 
at  ri  and  r2  for  an  emitter  ra  is 


F(ri,r2) 


/ 


Iagl92  ^ a, 


(3.1) 


where  gf  is  the  Green’s  function  for  the  propagation  from  emitter  a  to  the  receptor  i  (Cornwell, 
2004a,  equation  15). 


The  Green’s  function  for  the  propagation  of  the  radiation  from  emitter  a  to  receptor  1  is 
(Jackson,  1975,  Chapter  6) 

<Ji  =  t - — t  exp(ik\r a  -  ri  +  Dn\).  (3.2) 

|ra  -  ri  +  Dn\ 

Here  k  =  2ir/\  is  the  wavenumber. 

Expansion  of  the  product  of  the  Green’s  functions  (Cornwell,  2004a,  equation  18)  yields,  to 
hrst  order, 

V (ri,  r2)  =  j  /(ri,r2;  ra;  D)Iaexp[ik(rf  —  r%)/2D  —  ik(ji  -  r2)  •  ra/D\  dra  (3.3) 
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(cf.  Cornwell,  2004a,  equation  19).  Here  /(iq,  r2;  ra ;  D )  is  a  function  that  incorporates  the  amplitude 
factor  from  the  Green’s  functions. 

The  traditional  (u,v,w)  coordinate  system  of  radio  astronomy  has  u  =  k(ri  —  r2)/27r.  In 
standard,  far-held  radio  astronomical  imaging,  a  “direction  cosines”  coordinate  system  is  used  in 
which  tan  i  ~  sin£  ~  £  ~  xa/D  and  tanm  ~  sinm  ~  m  ~  ya/D.  We  relax  the  small-angle 
approximation  here,  but  to  emphasize  the  similarity  with  standard  radio  astronomical  imaging,  we 
use  the  coordinate  system  £  =  (£,m)  =  r a/D. 

Equation  3.3  can  thus  be  re-written  in  the  suggestive  form 

U(rq,r2)  =  j  flaex.p[ik(rl  —  r%)/2D  —  2iri(u£  +  vm)\dra 

=  j  Dfla  exp[ik(rf  -  rl) /2D\e~2^+vm)  d£dm,  (3.4) 

As  discussed  by  Carter  (1988),  a  redefinition  of  the  antenna  locations  results  in  the  first  exponential 
term,  the  Fresnel  term,  disappearing.  The  redefinition  applies  to  the  z  coordinate  of  the  antenna 
locations: 

z'  =  D-  y/D2  -  ( x2  +  y 2).  (3.5) 

Cornwell  (2004a)  calls  this  coordinate  redefinition  “Carter  refocussing.” 

Cornwell  (2004a)  continues  to  expand  the  phase  term  of  the  Green’s  function  product  (eqns.  [3.2] 
and  [3.3])  to  second  order  (see  Cornwell,  2004b).  After  “Carter  refocussing”  coordinate  redefinition, 
the  result  is 

V (ri,  r2)  =  ///.expHMn  “  r2)  '  ra/D]  exp [~ik(rj  -  rj)rl/2D3}  dra 
=  f  Dfla  exp[— ifc{(rf  -  rj) /2D}{ra/ D}2}e~2^ue+vm)  dldm. 

(3.6) 


Equation  (3.6)  clearly  has  the  form  of  a  2-dinrensional  Fourier  transform.  Further,  it  is  similar 
in  appearance  to  the  radio  astronomical  re-projection  expression  (Cornwell  et  al.,  2008,  eqns.  10-14). 
In  an  effort  to  make  the  similarity  even  more  apparent,  we  write  equation  (3.4)  as2 

V (ri,  r2)  =  J  N(n,r2-,D)(DfIa)e-2^ue+vmU£dm 

=  N(ri,r2-,D)<8>V(ri,r2\z1  =  0,Z2  =  0).  (3.7) 

where  <8>  is  the  convolution  operator. 

4.  DEPTH  OF  FOCUS 

Using  the  equivalent  of  the  re-projection  term,  in  the  N  function,  we  can  expand  the  phase 
term  to  obtain  an  estimate  of  the  depth  of  focus  or  the  range  in  distance  D  over  which  equa¬ 
tion  (3.7)  applies.  Equivalently,  this  “depth  of  focus”  also  implies  an  angular  distance  over  which 

2While  Cornwell  (2004a)  uses  G  to  denote  the  convolution  function,  we  prefer  N,  both  to  distinguish  it  from  the 
original  Green’s  functions  as  well  as  to  emphasize  the  near-field  aspects. 
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equation  (3.7)  applies.  Consider  an  emitter  at  position  ra  (Figure  1),  such  that  it  is  perpendicular 
to  Dh  (i.e.,  za  =  0).  This  emitter  can  be  considered  to  be  either  at  an  angular  offset,  with  a 
direction  cosine  given  by  ra/D  or  at  a  distance  \Jr2  +  I)2.  (The  latter  also  implies  a  different  h.) 

We  begin  with  just  the  phase  term 

‘^(S)2  =  -  ^VAD,  (4.1) 

where  q2  =  £2  +  m2  (Cornwell,  2004a,  eqn.  25).  The  quantity  \/r2  —  r2  can  be  no  larger  than  the 
maximum  baseline  in  the  array,  and  we  shall  approximate  r\  ~  rl  ~  B 2 .  Further,  the  Fresnel  zone 
is  defined  as  Rp  =  \J A D. 


We  define  the  depth  of  focus  as  the  change  in  the  distance  D  that  produces  a  phase  change  of 
order  1  radian.  Letting  D  — ►  D  +  AD  in  equation  (4.1),  expanding  to  lowest  order  in  A  D/D,  and 
setting  the  resulting  perturbation  term  equal  to  unity,  algebraic  manipulation  yields 


AD  ~ 


(4.2) 


Physically,  the  depth  of  focus  is  related  to  two  factors.  The  first  is  related  to  the  size  of  the  array 
relative  to  the  distance  to  the  target  D/B.  The  second  is  related  to  the  size  of  the  area  being 
imaged  relative  to  the  size  of  the  Fresnel  zone  Rpjra. 


5.  IMPLEMENTATION 


As  Cornwell  (2004a)  describes,  the  equivalent  convolution  function  can  be  written  as  a  Hankel 
transform.  In  order  to  make  this  explicit,  we  express  the  convolution  function  as 


N  =  J  ex.p[ik(rj  -  r%)q2 /2D}e2ni(ul+vmUldm 

=  [  exp [ik{r\  -  r2)g2/2L)]e2«96(cosxcoS^+sinxsm^)  dqd^ 


(5.1) 


Here  we  express  the  array  baselines  ( u ,  v )  and  the  target  direction  cosines  ( l ,  m)  in  a  polar  coordinate 
system,  (u,v)  =  (6  cos  i/j,  6  sin  ip)  and  ( l,m )  =  (g  cosy;,  9  sin  x).  Using  trignometric  identies,  we  can 
perform  one  of  the  integrations  (over  x),  leaving 


N  = 


exp [ik(r\  —  r2)q2 /2D]Jo(27rqb)  qdq, 


(5.2) 


where  Jq(x )  is  the  zeroth-order  Bessel  function.  The  left-hand  side  of  equation  (5.2)  is  now  in  the 
form  of  a  zeroth-order  Hankel  transform  of  the  (radially-symmetric)  function  exp[ifc(r^  —  r2)q2 /2D]. 
Implementations  of  the  Hankel  transform  are  available.3 


An  important  point,  though  often  left  implicit  in  Cornwell  (2004a),  is  that  the  visibilities  are 
expressed  in  terms  of  the  sky  brightness.  For  a  given  sky  brightness,  one  can  Fourier  transform 
(2-dinrension)  to  obtain  the  visibility  function,  then  re-project  it  to  form  the  full,  3-dimensional 

3e.g.,  http://wwwuser.gwdg.de/%7Emleuten/MATLABToolbox/HankelTransform.html 


6 


visibility  function.  This  process  is  limited  only  by  numerical  errors  (Cornwell  et  al.,  2008).  For 
imaging,  we  have  measured  the  visibilities  and  desire  the  (unknown)  sky  brightness. 

We  return  to  equation  (3.7).  Generally,  one  would  like  to  use  a  Fast  Fourier  transform  (FFT) 
to  obtain  the  image  I  from  the  visibility  V.  In  the  radio  astronomical  case,  the  visibility  data 
are  gridded  onto  a  uniform  grid  in  order  to  enable  the  FFT.  Consequently,  we  now  discuss  the 
implementation  in  matrix  form. 

V  (ri,  r2)  =  N(ri,r2;D)<8>V(r1,r2\zi  =  0,Z2  =  0) 

NTU({r;})  =  IVTIV-(8)y({ri}|{«i  =  0}).  (5.3) 

For  a  unitary  matrix,  NTN  =  I,  implying  that  the  r.h.s.  of  eqn.  (5.3)  reduces  to  just  F({r.;}|{zj  = 
0}). 


Implementation  of  near-field  rc-projection  thus  involves  the  following  steps,  for  a  given  array. 
We  assume  that  the  antenna  locations  r*  are  known  and  that  the  array  output  are  the  visibilities 
V (rt.  rf)  for  all  unique  pairs. 

1.  Choose  a  distance  at  which  the  target  is  known  or  suspected  to  be.  This  distance  should  be 
informed  by  the  depth  of  focus  (§4). 

2.  Compute  the  transpose  matrix  IVT.  This  matrix  is  independent  of  the  visibility  data,  so  that 
for  an  array  with  fixed  antenna  positions,  pre-computation  of  this  matrix  is  possible. 

3.  Grid  the  observed  three-dimensional  visibilities.  In  the  radio  astronomical  case,  this  is  per¬ 
formed  via  a  multiplication  of  the  observed  visibilities  by  a  gridding  convolution  function. 
The  gridding  convolution  function  is  chosen  to  have  a  small  support,  and  the  gridded  visi¬ 
bility  function  can  be  easily  resampled  on  a  uniform  grid.  Again,  it  may  be  possible  for  this 
function  to  be  pre-computed. 

4.  Multiply  the  observed  (and  gridded)  visibilities  by  the  IVT  matrix  to  produce  the  re-projected 
visibilities. 

5.  Perform  a  2-dinrensional  transform  of  the  near-field  re-projected  visibilities  to  obtain  the  sky 
brightness. 


As  Cornwell  (2004a)  and  Cornwell  et  al.  (2008)  discuss,  depending  upon  the  complexity  of  the  sky 
brightness,  additional  processing  steps  can  be  conducted  such  as  deconvolution  (e.g.,  clean). 

6.  APPLICATION  OF  NEAR-FIELD  W-PROJECTION  IMAGING  TO  THE  SPACE 
FENCE  (AFSSS) 

In  this  section  we  illustrate  the  use  of  near-field  re-projection  using  the  Space  Fence  receiving 
arrays.  Briefly  the  receiving  arrays  consist  of  St.  Andrew’s  crosses  with  a  maximum  east-west 
(E-W)  baseline  of  1600  ft.  (488  nr)  and  a  maximum  north-south  (N-S)  baseline  of  1200  ft.  (366  nr). 
The  observation  frequency  is  216.98  MHz,  producing  an  angular  resolution  in  the  E-W  direction  of 
3.4  nrilliradians  (11. 7')  and  in  the  N-S  direction  of  4.5  nrilliradians  (15. O').  Using  equation  (11)  from 
Carter  (1988),  the  far  field  is  at  least  1100  knr  distant  for  the  E-W  baselines  and  600  knr  distant  for 
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Fig.  2 — The  Long  Wavelength  Demonstrator  Array.  The  dipoles  are  in  the  foreground  and  are  approximately 
1  m  high.  The  LWDA  electronics  are  housed  in  the  shed  in  the  background.  The  operational  frequency  range 
is  60-80  MHz.  Antennas  of  the  the  Very  Large  Array  (VLA)  are  in  the  background. 


the  N-S  baselines.  For  additional  details  about  the  Space  Fence  transmitting  and  receiving  arrays, 
see  Meekins  (1988)  and  Solomon  (1990). 

We  begin  by  illustrating  the  concept  of  using  imaging  with  a  phased  dipole  array  to  localize 
emitters  within  the  field  of  view.  NRL  has  implemented  the  Long  Wavelength  Demonstrator  Array 
(LWDA),  an  array  of  16  dipoles  located  in  central  New  Mexico4  and  operating  over  the  frequency 
range  60-80  MHz  (Figure  2).  The  array  is  designed  as  a  testbed  for  radio  astronomical  observations, 
with  a  pseudo-random  distribution  of  dipoles  and  an  operational  frequency  range  chosen  to  be 
compatible  with  one  of  the  operational  frequencies  for  the  nearby  Very  Large  Array  (VLA). 

It  has  long  been  known  that  as  a  meteor  enters  the  Earth’s  atmosphere,  its  velocity  is  high 
enough  that  it  can  produce  a  ionized  trail,  which  can  reflect  a  transmitted  RF  signal  (e.g.,  Millman 
et  al.,  1948).  In  2006  November,  the  LWDA  acquired  data  at  a  frequency  of  61  MHz  during  the 
Leonid  meteor  shower.  This  frequency  was  chosen  to  be  one  at  which  a  number  of  TV  stations 
in  New  Mexico  and  the  southwest  U.S.  transmit.  Figure  3  illustrates  two  examples  of  a  meteor 
reflection  imaged  by  the  LWDA. 

LWDA  imaging  of  meteor  trail  reflections  bears  obvious  similarities  to  the  Space  Fence  receiv¬ 
ing  sites.  Both  are  phased  dipole  arrays,  with  operational  frequencies  differing  by  a  factor  of  only 
3.6,  and  both  can  image  the  reflection  of  a  signal  from  a  distant  transmitter.  Because  it  is  intended 
to  test  radio  astronomical  hardware  and  software,  the  LWDA  images  are  formed  by  Fourier  trans¬ 
forming  the  visibility  data  (phase  differences)  from  unique  dipole  pairs.  Even  for  reflections  from 
meteor  trails,  this  far-held  imaging  approach  is  sufficient  as  the  far-held  distance  is  at  an  altitude 
of  merely  0.5  km.  In  contrast,  the  higher  frequency  and  larger  size  of  the  Space  Fence  receiving 

4Latitude:  34.076N,  Longitude  107.618W;  west  of  the  city  of  Socorro,  NM 


Table  1 — Space  Fence  Depth  of  w-projection  Focus  (E-W) 

300  km  altitude  1000  km  altitude 


Transverse 

distance 

(km) 

Angle 

(rad.) 

Field 

of  View 

(°) 

AD 

(km) 

Field 

Angle 

(rad.) 

of  View 

(°) 

AD 

(km) 

10 

0.03 

3 

300 

30 

0.10 

11 

300 

0.03 

3 

1000 

100 

0.33 

38 

300 

0.10 

11 

1000 

300 

1.00 

115 

300 

0.30 

38 

1000 

445 

1.48 

170 

238 

1000 

1.00 

115 

1000 

1483 

1.48 

170 

1000 

arrays  means  that  the  far-held  distance  is  at  a  higher  altitude  and  near-held  effects  must  be  taken 
into  account  in  making  the  images. 

Table  1  summarizes  the  depth  of  focus  (equation  4.2)  for  the  Space  Fence  receiving  (east-west) 
arrays  for  satellites  at  two  different  altitudes.  Several  comments  are  warranted  about  this  table. 
First,  many  of  the  depth  of  focus  values  AD  are  equal  to  the  assumed  altitude  D.  This  equality 
is  the  result  of  our  assumption  in  deriving  eqn.  (4.2)  that  AD  is  “small”  so  that  an  expansion  in 
terms  of  A  D/D  is  valid.  Allowing  A  D/D  — >  1  may  seem  to  violate  this  assumption,  but  we  have 
made  other  approximations  (e.g.  ,  r\  -  r\  ~  B 2)  that  are  conservative. 

Second,  examination  of  Figure  1  shows  that  the  depth  of  focus  applies  to  an  angular  displace¬ 
ment  from  the  center  of  the  held  ( ra/D ).  Consequently,  the  full  held  of  view  (Columns  3  and  6)  or 
the  area  that  can  be  imaged  with  one  depth  of  focus  correction  is  twice  the  angular  displacement. 
Further,  anticipating  the  use  of  the  Space  Fence  to  large  helds  of  view,  we  extend  the  list  of  trans¬ 
verse  displacements  to  correspond  to  an  angular  displacement  of  1.48  radians  (85°  oh  zenith,  or  a 
held  of  view  of  170°). 

Our  interpretation  of  these  values  is  the  following.  At  altitudes  of  about  1000  km,  the  depth  of 
focus  for  rc-projection  imaging  is  comparable  to  the  altitude,  even  for  large  helds  of  view.  Thus,  a 
single  ir-projection  should  be  able  to  image  most  of  the  near-held  region.  The  one  exception  would 
be  for  low-altitude  satellites  (^  300  km)  at  large  angles  oh  zenith.  For  these  directions,  one  may 
need  to  refocus  the  array  (i.e. ,  choose  a  different  direction  for  n,  which  in  turn  implies  a  different 
set  of  [ u ,  v ,  w]). 

A  nominal  CONOPS  for  near-held  w- projection  imaging  with  the  Space  Fence  is  thus  the 
following.  Phase  data  are  acquired  as  in  current  standard  operating  procedure.  Visibility  data5 
(i.e.,  correlation  coefficients)  are  constructed  by  determining  the  phase  differences  between  the 
unique  pairs  of  antennas,  relative  to  a  (arbitrary)  reference  antenna.  The  visibility  data  are  then 
processed  in  two  separate  methods  or  “streams”  (which  could  clearly  be  conducted  in  parallel): 

'In  general,  visibility  data  are  complex,  having  amplitude  and  phase  or  real  and  imaginary  components.  In 
practice,  for  the  Space  Fence,  only  the  phase  data  are  measured,  and  it  is  sufficient  to  choose  a  nominal,  constant 
amplitude  for  all  antennas. 


Near-field  w-projection:  NRL/MR/7210— 09-9173 


9 


Near-field  The  prescription  of  Section  5  is  used,  e.g.,  with  a  target  altitude  of  1000  km.  In 
particular,  because  one  depth  of  focus  is  sufficient  and  the  positions  of  the  Space  Fence 
receiving  array  elements  are  fixed,  Steps  #1  and  # 2  can  be  precomputed  (e.g.,  stored  in  a 
look-up  table).  Further,  the  gridding  function  of  Step  #3  could  also  be  precomputed. 

Far-field  The  visibility  data  are  gridded  and  transformed  (FFT).  Again,  because  the  positions  of 
the  receiving  arrays  are  fixed,  the  gridding  function  could  be  precomputed. 


The  output  of  either  processing  stream  would  be  an  image,  showing  energy  as  a  function  of  position 
above  the  array  (Figure  3).  The  image  is  effectively  a  2-dimensional  matrix,  the  elements  of  which 
encode  the  direction  cosine  information  to  any  radar  reflections  within  the  field  of  view.  Quality 
control  checks  can  be  performed  on  such  an  image.6  An  image  that  passes  the  quality  control 
check(s),  could  then  be  scanned  in  an  automatic  manner  for  any  energy  intensity  peaks  representing 
reflections  from  satellites  passing  through  the  Space  Fence  beam.  The  position  of  the  intensity  peaks 
would  then  be  extracted  from  the  image  and  used  to  report  satellite  positions. 

Finally,  we  note  a  curious  fact  about  near-field  processing  for  the  Space  Fence.  One  of  the 
original  motivations  for  the  w-projection  technique  was  that  radio  astronomical  interferometers  are 
not,  in  general,  coplanar,  because  of  local  topology,  the  array  is  large  enough  that  the  radius  of 
the  Earth  has  to  be  taken  into  account,  or  both.  The  Space  Fence  is  coplanar  (Solomon  1990; 
E.  Lydick  2009,  private  communication),  but,  in  order  to  use  it  for  near-field  imaging,  it  must  be 
transformed  into  a  non-coplanar  array  (viz.  the  discussion  concerning  eqn.  [3.5]). 


6For  example,  a  standard  radio  astronomical  test  is  to  measure  the  value  of  the  rms  pixel  intensity  in  the  image. 
Only  images  with  rms  pixel  intensity  values  lower  than  some  threshold  are  considered  viable  for  further  processing. 
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Fig.  3 — LWDA  images  of  reflected  TV  signals  from  ionized  meteor  trails  during  the  2006  Leonid  meteor 
shower.  Each  image  is  constructed  from  a  0.1  s  data  acquisition  at  a  frequency  of  61  MHz.  ( Top  Left) 
and  ( Top  Right)  Images  of  reflected  signals  from  Leonid  meteors.  ( Bottom )  A  reference  image  showing  the 
nominal  sky  observed  by  the  LWDA.  The  zenith  is  at  the  center  of  the  image.  The  two  discrete  emitters 
near  the  center  of  the  field  are  celestial  radio  sources  (Cas  A  and  Cyg  A) .  The  band  of  emission  stretching 
from  the  upper  right  to  lower  left,  terminating  in  the  strong  emitter  near  the  edge  of  the  field  is  from  the 
Milky  Way  Galaxy. 
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